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Abstract: We show that a cavity optomechanical system formed by a 

mechanical resonator simultaneously coupled to two modes of an optical 
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scheme is realizable for sufficiently strong single-photon optomechanical 
coupling in the resolved sideband regime, and is robust against cavity losses. 
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1. Introduction 

Simple quantum information protocols and quantum gates have been recently implemented 
with high fidelity in trapped ions and in circuit cavity QED (see e.g.. Ref. Q] for a review). In 
an efficient quantum network, the information elaborated by a solid state processor at a node 
should be then robustly encoded in single-photon qubits for long-distance communication and 
distribution of quantum information. The possibility to implement high-fidelity two-qubit gates 



between single-photons would greatly facilitate such quantum information routing; an example 
is provided by perfect Bell-state discrimination for quantum teleportation and entanglement 
swapping 0, which could be implemented deterministically if a quantum phase gate (QPG) 
between single photon qubits would be available 0 . It is well known that to obtain such a QPG 
one needs nonlinearities implicit in post-processing measurements or a nonlinear system 0. 
A particularly convenient QPG is obtained when the conditional phase shift between the two 
photonic qubit is equal to n, because in this latter case the QPG is equivalent, up to local 
unitary transformations, to a C-NOT gate 02). Using all-optical solutions this is not easy to 
achieve because to process the information one needs strong photon-photon interaction. In fact, 
to implement quantum information with photons, a nonlinear interaction is needed either to 
build a two-photon gate operation © (but the commonly achievable y? susceptibility factor is 
typically too small), or the nonlinearity implicit at the detection stage in linear optics quantum 
computation 0. 

Then, one has to think differently as for example using electromagnetic ally induced trans¬ 
parency (EIT) |U[9). Indeed the possibility of reducing the speed of light in atomic media 
with A-type levels was proposed and experimentally obtained ED- Although the two weak 
fields wave functions traveling with slow group velocity can have a very high nonlinear cou¬ 
pling when they propagate in an atomic media with A-like level configurations HTHTD . we 
need, however, to obtain the same result at the single photon level and this is extremely hard 
to achieve or even impossible with A-like configurations. A full quantum analysis has shown 
that a large nonlinear cross-phase shift is achievable using an atomic M level structure (131415) . 
but it was also shown that a trade-off between the size of the conditional phase shift and the 
fidelity of the gate exists. This can be avoided in the transient regime, which is however exper¬ 
imentally challenging. More recently, important results have been achieved by exploiting two 
different solutions able to provide the required effective nonlinearities: i) the strong dispersive 
coupling of light to strongly interacting atoms in highly excited Rydberg states HMD; ii) a 
fiber-integrated cavity QED system employing a whispering gallery mode resonator strongly 
coupled to a single Rubidium atom (22) . 

Nevertheless, in recent years it has been proposed lf23l . and then experimentally real¬ 
ized lf24l427l . that EIT-like effects could be obtained also within cavity optomechanical sys¬ 
tems. Furthermore, it is well known that the ponderomotive action of light, together with the 
backaction of the mechanical oscillator interacting with it, is responsible for an effective optical 
Kerr nonlinearity (28), which in turn, may give rise to interesting quantum phenomena, such as 
squeezing of the cavity output light, as predicted a couple of decades ago 112911301 and recently 
experimentally achieved (311433) . 

In this paper we will show that a QPG for simple photonic qubits with a conditional phase 
shift equal to n is achievable by employing a cavity optomechanical system with sufficiently 
large single-photon optomechanical coupling, and relatively high mechanical Q-factor Q m . 
Refs. (341135 1 first suggested that multi-mode optomechanical systems in the single-photon 
strong coupling regime could be exploited for quantum information processing with photons 
and phonons, and a first example of optomechanical implementation of a QPG has been re¬ 
cently provided in Ref. (36) . Here we further develop these ideas, by proposing a much simpler 
scheme, which requires the control of only two cavity modes and of a single mechanical res¬ 
onator, rather than four optical cavity modes and two mechanical resonators as in Ref. (36) . 
Recent progress in the realization of strongly coupled nano-optomechanical systems B3H2D 
suggests that the QPG scheme proposed here could be implemented in the near future. 


2. The Model 


We consider an optomechanical system consisting of a mechanical resonator interacting with 
two optical modes, which is described by the following Hamiltonian 

H = TiO)\a\di+Tio>2.a\a2 + Tico m b'b + Ti(^g\a\a\-\-g2d\a2 S j (b + P), (1) 

where, dj(b) and a\ (b 1 ) are the annihilation and creation operators for the optical cav¬ 
ity (mechanical) modes, with frequency a>i/2n and co m /2n respectively, and with \d,,d] = 
[i>,P] = 1 ; gi = (dcOi/dx)x z p f is the i-th single-photon optomechanical coupling rate, with 
x zp f = y/fi/2m(O m the spatial size of the zero-point fluctuation of the mechanical oscillator. 

We focus on the simplest choice for an optical qubit, the state space spanned by the lowest 
Fock states of an optical mode, |0) and 11); to be more specific we want to implement a QPG 
between the optical qubits associated with two optical cavity modes of the optomechanical 
system under study. The generic initial (pure) state of the two optical qubits is given by 

I w)m = o^o|0)i|0)2 + «oi|0)i|1)2 + o:io|1)i|0)2 + o:ii|1)i|1)2j (2) 

corresponding in general to an entangled state of the two modes with up to two photons. A 
simple proof-of-principle demonstration could be achieved by restricting to factorized input 
states of the two cavity modes, which could be provided by two weak laser pulses driving the 
two selected cavity modes at frequencies (On and ft)/ 2 , similar to the preliminary experimen¬ 
tal demonstration of a QPG given in Ref. |5). In this case the two cavity modes are prepared 
in a product of two coherent states with amplitudes a c n and Gfe, 1^(0)) = |®rfi)i I®tf 2)2 — 
[|0)i +a f /i|l)i] [|0)2 + C^/2|1)2], where the latter expression is valid for |a^i|, |« z / 2 1 -C 1. For- 
input laser powers Pj, cavity detunings A j = C0j — and decay rates Kj, j = 1,2, the ampli¬ 
tudes are given by a ( ij = 

It is convenient to move to a frame rotating at the corresponding driving laser frequency 
for each cavity mode, providing therefore the phase reference for each optical qubit; this is 
equivalent to move to the interaction picture with respect to the free optical Hamiltonian Hq = 
h(Oi\d[d\ J r h(Oi 2 d\d 2 , in which the system Hamiltonian becomes 

H = Mini +TiA2ti2 + 1i(O m b , b + Ji(o m fn (b + P), (3) 

where we have used the cavity mode photon number operators fij = d-dj, and we have defined 

fn = [g\d\ +g2«2]M„. 
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3. Hamiltonian dynamics 

In order to have a physical description of how the effective optical nonlinearity provided by the 
optomechanical interaction allows to implement the QPG, we first study the ideal case with no 
optical and mechanical losses, in which the dynamics is determined solely by the Hamiltonian 
of Eq. d3}, i.e., by the unitary operator 0 (t ) = e ,//f +. In such a case the dynamics can be ex¬ 
actly solved: in fact, profiting from the fact that both photon number operators fij are conserved, 
and moving to a photon-number-dependent displaced frame for the mechanical resonator, one 
can rewrite the unitary evolution operator in a form in which the optical and mechanical evo¬ 
lution operators are conveniently factorized. In fact, acting with the photon number conserving 
mechanical displacement operator 

D(f„) = exp {(b 1 -b)ffi\ , 


(4) 






which separates the Hamiltonian according to 


o (fn) HD' (fit) =H opt + H b , 


(5) 
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one gets 

0(0 = ^opt(0O + (/*) 0 6 (O0 (/*), (9) 

where U opt (t) = and U b (t) = e ~ lHbt ' h . Therefore the dynamics of the two optical 

modes is mostly determined by the effective unitary operator U opt (t), possessing either self-Kerr 
terms « rij and the cross-Kerr term —2ftgig 2 tnin 2 /ft) m ; the other factor D 1 (/«) O b (t)D (fn) 
however also affects the optical mode dynamics since it entangles them with the mechanical 
resonator, by correlating the resonator position with the photon numbers. 

The photon number conserving dynamics allows to stay within the logical space described 
above, i.e., the one spanned by optical Fock states with no more than one photon (and this 
will remain true even when we will include optical losses). In this case one can always fix 
the two detunings in order to eliminate completely the effect of the self-Kerr terms. In fact, 
within this subspace, nj = hj, and therefore, taking Ay = g 2 -/co m , j = 1,2, one has the effective 
unitary operator U op t(t) = exp[2igig2tnin2/(O ni ], yielding a nonlinear conditional phase shift 
0ni (t) = 2 .g\g 2 t/(O m only when each cavity mode has one photon, i.e., 


2sIS?' 
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( 10 ) 

Therefore we expect to get a conditional phase shift equal to n when the interaction time t is 
equal to 


71 COjji 
2gigi' 


( 11 ) 


We now evaluate the exact Hamiltonian evolution in order to see to what extent the in¬ 
teraction with the mechanical resonator affects and modifies the ideal QPG dynamics de¬ 
fined by Eqs. Cl and CD- The natural choice for the initial state is the factorized state 
p (0) = | y/)m(ys\ <S> p*, where | (//),■„ is the generic initial state of Eq. (H, and p* is the thermal 
equilibrium state of the mechanical resonator, with h mean thermal phonons. By renumbering 
|0) 1 10>2 -> |0>, 11) 1 10) 2 -f 11), |0)i|l) 2 -4- \2), 11) 1 11) 2 -4- |3), we can write the state of the 
whole system at time t as 


p(t) = U(t)\y) in (v\®ptuHt) ( 12 ) 

= £ cc k a?U opt (t)\k) (A) U b {t)D (f k ) p*D + (/,) ^ (t)D (f,), 
k,l =0 


where 

D(f k )=exp[f k (b r -b)] (13) 

is now a displacement operator acting only on the mechanical resonator degree of freedom, 
with a c-number displacement f k , with f Q = 0, f\ = gi/(O m , f 2 = gi/(Om, h = (gi + g 2 )/®,„. 





We are interested in the state of the two optical qubits only, and therefore we have to trace over 
the mechanical resonator. Using the explicit expression of t/ op t(f) (with the choice of detuning 
specified above), and performing the trace, the reduced state of the optical modes reads 
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where <5* j is the Kronecker delta, while Ck.i(t) is the factor describing the decoherence caused 
by the interaction with the mechanical resonator and whose explicit expression is given by (see 
the appendix for its derivation) 


Ck,i{t) =exp 


(fk - fif (1 - cos (O m t) (2n + 1) + i (f? - fl ) sin (o m t 


(15) 


We quantify the QPG performance with the fidelity relative to the ideal pure target state corre¬ 
sponding to a n conditional phase shift, that is 



(16) 


The corresponding fidelity F(t) can be written as 
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Actually, the QPG performance can be characterized by the gate fidelity, which is the average 
of the above quantity over all possible input states of the two qubits ED- Since the averages 
are given by |a A | 4 = 1/8, \/k, and |a*| 2 |a/| 2 = 1 /24, \/k fi l ED . using the explicit expression 
for c* /(f) of Eq. (ITD implying in particular that c* A (f) = 1 Vk and that c*/(f) = c/. A (f)*, and 
the explicit values of the c-numbers fk, we get 
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From Eqs. (TiTl i one can see that the gate fidelity achieves the ideal value of unity when 
two conditions are satisfied: i) the conditional phase shift is equal to 71, (or more generally 




































to an odd multiple of it), 2g\g2t/(O m = (2m + 1 )n (integer m); ii) Q./(f) = 1, \/k,l. Eq. (flSl ) 
shows that the latter conditions are achieved for generic nonzero couplings gi and g 2 only 
after every mechanical oscillation period, i.e., when (O m t = 2pJt, p = 1,2,.... The QPG will be 
minimally affected by losses for the shortest interaction time t„ of Eq. (fTTb . and therefore, the 
ideal conditions for a unit-fidelity QPG with a conditional phase shift equal to n are 

2 g 1 k-G ^ + Ott- /ioa 

— It (Ofntji — 2 K =$• g\g2 — . (19) 

CO,,, 4 

Therefore, when optical and mechanical losses are negligible, one can realize an ideal QPG with 
a simple optomechanical setup by fixing the interaction time between the mechanical resonator 
and the cavity modes according to Eq. dTTt . and provided that the single-photon optomechanical 
couplings can be tuned to the strong coupling condition of Eq. G3. The interaction time can 
be controlled in tunable optomechanical systems in which the interaction can be turned on and 
off, as it could be done for example in the optomechanical setup of Ref. [40|, where a vibrating 
nanobeam is coupled to the evanescent field of a whispering gallery mode of a microdisk. It 
is also relevant to stress that under these conditions the QPG is practically insensitive to the 
effect of thermal noise acting on the resonator, because at the exact gate duration t K the fidelity 
becomes completely independent upon the mean thermal phonon number n. Eq. (IT8l) shows that 
instead the gate fidelity significantly drops for increasing n as soon as t ^ t K . 

4. Dissipative dynamics 

Let us now consider the realistic situation in order to see to what extent optical losses, mechan¬ 
ical damping, and thermal noise affect this ideal gate behavior. In that case the evolution is no 
more analytically tractable, and we will consider the numerical solution of the master equation 
for the density matrix of the optomechanical system under study. 

Introducing the cavity modes decay jsj (i = 1,2), the mechanical damping y,„ = COm/Qm , and 
the mean thermal phonon number associated with the reservoir of the mechanical resonator, n, 
the master equation in the usual Born-Markov approximation can be written as ll42l 

j t P (0 = ^fi,p(t)\ + ^{ 2 aip(t)d\-a\dip(t)-p(t)a\ai) 

+ -y ( 2 d 2 p( 0«2 “ al& 2 p(t) - p(t)a\& 2 ) (20) 

+-^p(n + 1)(2 bp(t)b' — b'bp(t) — p(t)b'b) + ^-n(2b^ p(t)b — bP p(t) — p(t)bb'), 

where H is the Hamiltonian in Eq. ({3}. 

In Fig. 1 we compare the behavior of the gate fidelity J^(f) in the absence of damping and 
losses of Eq. (fl8l) versus the dimensionless interaction time (o m t , either at n = 0 (red dot-dashed 
curve), and at n = 10 (full black curve) with the corresponding curves in the presence of optical 
and mechanical damping processes. These latter curves (the dotted blue line corresponds to 
n = 0, while the green dashed line to h= 10) are obtained from the numerical solution of the 
master equation Eq. ( l20t in the case K‘| = Ki = 10 2 0 ) m .Q, n = 10 6 , while we have fixed the 
couplings according to the ideal strong coupling condition of Eq. ( IT9b g\ = g 2 — ft),„/2. 

We see, as expected, that in the absence of optical and mechanical losses, &(t) = 1 exactly 
at the interaction time t K of Eq. CD, regardless the value of the temperature of the mechanical 
reservoir. At different times the gate performance is strongly affected by thermal noise; what 
is relevant is that in the presence of realistic values of mechanical damping and of optical loss 
rates, this scenario is still maintained, with a limited decrease of the gate fidelity. 




Fig. 1. Numerical solution of the master equation Eq. 120b for the gate fidelity &{i) 
versus the dimensionless interaction time co m t. We compare four different cases: i) zero 
damping and losses y m — K\ = K 2 = 0 and n = 10 (full black line); ii) zero damp¬ 
ing and losses and n = 0 (red dashed-dotted line); iii) with damping and losses (JCj = 
Ki = 10 “U) m , Q,„ = 10 6 ,) and n = 10 (green dashed line); iv) with damping and losses 
()Ci = K2 = IQ - 2 CO,,,. Q m = 10 6 ,) and ii = 0 (blue dotted line). The numerical solutions 
for the zero damping and loss case are indistinguishable from the analytical expression of 
Eq. i(T 8} either at n = 0 and at n = 10. In all cases we have fixed the couplings according to 
the ideal strong coupling condition of Eq. {UK gl=g2 = C0m/2. 


5. Conclusions 

We have proposed a simple optomechanical setup which is able to implement an ideal QPG 
with a conditional phase shift equal to K between two optical qubits associated with the lowest 
Fock states (zero and one photon) of an optical cavity mode. The scheme is minimal because it 
employs only two modes of a high-finesse optical cavity and a single mechanical resonator cou¬ 
pled to them. The scheme is robust in the presence of realistic values of optical and mechanical 
losses and, if the interaction time is appropriately fixed, it is almost completely insensitive to 
the thermal noise acting on the mechanical resonator. The most stringent and challenging con¬ 
dition is the required strong optomechanical coupling condition, given by Eq. < IT9b . in which 
the single-photon optomechanical coupling must be of the order of the mechanical resonance 
frequency. Such a condition has not been achieved yet in current solid-state nanomechani¬ 
cal setups, for which record values corresponds to gj/(O m ~ 10 " 3 lf37H39| . On the contrary, 
such a strong coupling situation is normally achieved in ultracold atom realizations of cav¬ 
ity optomechanics, where the mechanical resonator corresponds to the collective motion of an 
ensemble of trapped ultracold atoms; for example one has g/ co m ~ 0.3 in Ref. 071 . The limita¬ 
tion in these latter systems is represented by cavity losses, because in this case one is typically 
far from the resolved sideband regime k /( 0 „, <C 1 which is required here in order that cavity 
losses do not alter significantly the effective cross-Kerr nonlinear interaction mediated by the 
resonator, responsible for the QPG dynamics. Therefore the present proposal could be imple¬ 
mented in experimental optomechanical platforms able to combine a significantly large single 











photon coupling g/c0 m ~ 0.5 with a resolved sideband operation condition k/(O m <C 1. 


6. Appendix 

We now derive the explicit expression for the decoherence coefficients c k j(t) of Eq. ( 1 1 5| >. From 
Eqs. <fj~2~b - <n~4l > and using the cyclic property of the trace, one has 


Ck,l(t) = Tr, (f,)U^t)D(f,)D^ C f k )U b (t)D{f k )Pb 


( 21 ) 


which is a thermal average of a combination of displacement operators. We notice then that the 
factor Ufi(t)D ( fi)D T (/*) U b (t) within the trace is just the Heisenberg time evolution for a time 
t of a displacement operator of a free mechanical resonator, so that 

ut(t)D(f,)fr (f k )U b (t) = exp [(Pe^-be-^) (/,-/*)] . (22) 

Inserting this solution within Eq. (|2TT) . and using the property of the displacement operator 
D(a)D(p) =D(a + /3)exp[iIm(a/3*)], one gets 


Ck,i(t) = Tr b {D - 1)] p*} exp [i{ff -fj;)sm(o m t\ . 

Performing the thermal average of the final displacement operator, according to 


(exp [ab' 


a*b] )th = exp 



we finally get Eq. ( fl~5l >. 


(23) 


(24) 
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